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Abstract 

We consider the general Landau-Lifshitz-Gilbert (LLG) dynamical theory underlying the mag- 
netization switching rates of a thin film uniaxial magnet subject to spin-torque effects and thermal 
fluctuations (thermal noise). After discussing the various dynamical regimes governing the switch- 
ing phenomena, we present analytical results for the mean switching time behavior. Our approach, 
based on explicitly solving the first passage time problem, allows for a straightforward analysis of 
the thermally assisted, low spin-torque, switching asymptotics of thin film magnets. To verify our 
theory, we have developed an efficient GPU-based micromagnetic code to simulate the stochastic 
LLG dynamics out to millisecond timescales. We explore the effects of geometrical tilts between 
the spin-current and uniaxial anisotropy axes on the thermally assisted dynamics. We find that 
even in the absence of axial symmetry, the switching times can be functionally described in a form 
virtually identical to the collinear case. We further verify that asymptotic behavior is reached fairly 
slowly, thus quantifying the role of thermal noise in the crossover regime linking deterministic to 
thermally assisted magnetization reversal. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 



More than a decade has passed since spin-torque effects were demonstrated experimen- 
tally by the switching of the magnetization of a thin ferromagnetic film when current is 
passed between it and a pinned ferromagnetic layer-"-. A spin-polarized current passing 
through a small magnetic conductor will deposit spin-angular momentum into the magnetic 
system. This in turn causes the magnetic moment to precess and in some cases even switch 
direction. This has led to sweeping advances in the field of spintronics through the devel- 
opment and study of spin- valves and magnetic tunnel junctions (see, for example,-). The 
theoretical approach to such a problem has conventionally been to treat the thin ferromag- 
netic film as a single macrospin in the spirit of Brown^. Spin-torque effects are taken into 
account phenomenologically by modifying the macrospin's LLG dynamical equation^. A 
thorough understanding of the phenomena, however, cannot proceed without taking into 
account the effect of thermal fluctuations. This is of particular experimental relevance since 
spin-transfer effects on nanomagnets are often conducted at low currents, where noise is 
expected to dominate. Recent debate in the literature over the proper exponential scaling 
behavior between mean switching time and current shows how the thermally assisted prop- 
erties of even the simplest magnetic setups leave much to be understood^"—. The interplay 
between spin-torque and thermal effects determine the dynamical properties of recent exper- 
imental studies on nanopillar devices^. Except at very high currents where the dynamics 
are predominantly deterministic, the switching appears to be thermal in nature. Fitting to 
experimental data requires accurate knowledge of the energetics which, in the realm of spin- 
torque, are hard to come by due to the inherently non-conservative nature of the spin-torque 
term. 

Theoretical progress has been hindered by the computational power needed to run numer- 
ical simulations to the desired degree of accuracy. The LLG equation, modified into its set 
of coupled stochastic equations can be studied in one of two ways: either by concentrating 
on the associated Focker-Planck equation or by constructing a stochastic Langevin integra- 
tor to be used enough times to gather sufficient statistics on the phenomenal. The latter 
approach, however, has been unable to extrapolate to long enough times to capture the dy- 
namical extent of the thermal regime. A recent paper by Taniguchi and Imamurai^, suggests 
that previous analytics of the thermally assisted dynamics should be revisited. Nonethe- 
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less, no numerical simulation has yet been able to evaluate the accuracy of the Taniguchi 
and Imamura results^. In our paper, we will show that simulations run harnessing the 
vast computational parallelization capabilities intrinsic in Graphics Processing Unit (GPU) 
technology for numerical modeling can allow a deeper probing of such a thermally activated 
regime. 

II. GENERAL FORMALISM 

A simple model of a ferromagnet uses a Stoner-Wohlfarth monodomain magnetic body 
with magnetization M. The body is assumed to have a size Im along the e^, direction, and 
size a in both the and directions. The total volume of the object is then V = a?lm- The 
shape of the body is assumed to be isotropic and the energy landscape experienced by M is 
generally described by three terms: an applied field H, a uniaxial anisotropy energy Uk with 
easy axis along yvk in the — making an angle u with the axis, and an easy-plane 
anisotropy Up in the Gy — plane, with normal direction fi/) perpendicular to yyk- The 
magnetization M is assumed to be constant in magnitude and for simplicity normalized into 
a unit direction vector m = M/|M|. A spin-polarized current J enters the magnetic body in 
the —Gy direction, with spin polarization factor 77, and spin direction along e^. The current 
exits in the same direction, but with its average spin direction aligned to that of M. The 
self-induced magnetic field of the current is ignored here as the dimension a is considered to 
be smaller than lOOOA, where the spin-current effects are expected to become dominant over 
the current induced magnetic field. The standard Landau-Lifshitz-Gilbert (LLG) equation 
used to describe the dynamics is then written as: 



where 7' = 7/(1+0^) is the Gilbert ratio, 7 is the usual gyromagnetic ratio and j = {h/2e)'r]J 
is the spin-angular momentum deposited per unit time with 77 = {J^ — Ji)/{J^ + J^) the 
spin-polarization factor of incident current J. The last two terms describe a vector torque 
generated by current polarized in the direction rip. These are obtained by assuming that 
the macrospin absorbs angular momentum from the spin-polarized current only in the 



m = — 7'm X Hcfj — a'-f'm x (m x Hcfj) 




(1) 
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direction perpendicular to m.- 

To write Hcs explicitly, we must construct a proper energy landscape for the magnetic 
body. There are three main components that need to be considered: a uniaxial anisotropy 
energy Uk, an easy-plane anisotropy Up and a pure external field interaction Uh- These are 
written as follows: 



Uk = -K{nK ■ 
Up = Kp{tvd ■ m)2 
Uh = -MsVm ■ H^^i 

In these equations, Ms is the saturation magnetization, Kp is the easy-plane anisotropy, 
K = {1/2)MsVHk and Hk is the Stoner-Wohlfarth switching field (in units of Teslas). In 
what follows, we consider the simplified model where we ignore the effects of easy-plane 
anisotropys and all external magnetic fields are absent. The full energy landscape then 
becomes f/(m) = Uk + Uh and reads: 



f/(m) = -K Uuk ■ m)2 + 2h ■ ml , (2) 

where h = Hext/HK, iiK is the unit vector pointing in the orientation of the uniaxial 
anisotropy axis, and effects due to external magnetic fields are included for reasons which 
will be explained in the next section. Such an energy landscape generally selects stable 
magnetic configurations parallel and anti-parallel to uk- The effective interaction field Hgff 
is then given by 

Hcff = -— ^V„,f/(m) = Hk [(fix ■ m)n^ + h] . (3) 

The symmetries of the problem lead to slightly simplified equations and the deterministic 
LLG dynamics can then be expressed as: 



rh = — m X [{n-K ■ m)nA' + h] — am x [m x {{Uk ■ m)nA' + h)] 
—aim X f m X k] + a^/m x k, 



(4) 
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where we have defined for future convenience / = j/{aHK) and introduced the natural 
timescale r = yHxt. 

III. THERMAL EFFECTS 

Thermal effects are included by considering uncorrelated fluctuations in the effective 
interaction fleld: Hcff — Hcff + Hth- These transform the LLG equation into its Langevin 
form. We model the stochastic contribution Hj/j by specifying its correlation properties, 
namely: 



(Rth) = 

{Htf,4t)HtHAt')) = 2D6,,kS{t - t') (5) 

The effect of the random torque Hth is to produce a diffusive random walk on the surface 
of the M-sphere. An associated Focker-Planck equation describing such dynamics was con- 
structed by Brown^. At long times, the system attains thermal equilibrium and, by the 
fluctuation-dissipation theorem we have: 

~ 2/i(l + a2) - 2(1 + ^ > 

It is convenient to introduce the notation K/ksT = C,, representing the natural barrier 
height of the uniaxial anisotropy energy. 

Setting now the external magnetic fleld to zero and considering only thermal fluctuations, 
the stochastic LLG equation reads: 



rfii = Ai{m) + Bikim) o Hth,k (7) 



where 



A(m) = al am x k — m x (^m x 

— {iiK ■ m) [m X hx — Oi {jik — ip-K ■ m)m)] 
Bik{m) = -eijknij - a{mimk - Sik). 
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and ^oHth,k means to interpret our stochastic dynamics in the sense of Stratonovich^I cal- 
culus in treating the multiplicative noise terms^^. 

We numerically solve the above Langevin equations by using a standard second order 
Heun scheme to ensure proper convergence to the Stratonovich calculus. At each time step, 
the strength of the random kicks is given by the fluctuation-dissipation theorem. Statistics 
were gathered from an ensemble of 5000 events with a natural integration stepsize of 0.01. For 
concreteness, we set the Landau damping constant a = 0.04. Different barrier heights were 
explored although the main results in this paper are shown for a barrier height of ^ = 80. To 
explore the simulations for long time regimes, the necessary events were simulated in parallel 
on an NVidia Tesla C2050 graphics card. To generate the large number of necessary random 
numbers, we chose a proven combination^^ of the three-component combined Tausworthe 
"taus88"— and the 32-bit "Quick and Dirty" LCG^. The hybrid generator provides an 
overall period of around 2^^^. 

IV. SWITCHING DYNAMICS 

In experiments, one is generally interested in understanding how the interplay between 
thermal noise and spin torque effects switch an initial magnetic orientation from parallel 
to anti-parallel and vice versa. The role of spin-torque can be clarified by considering how 
energy is pumped in the system, from an energy landscape point of view. As in the previous 
section, the magnetic energy of the monodomain is: 



The change in energy over time can be obtained after some straightforward algebra and is 
found to be: 



This expression shows how current pumps energy into the system. In the absence of current, 
the damping dissipates energy and, as one would expect, the dynamical flow is toward the 
minimum energy configuration. The sign preceeding the current term allows the expression 
to become positive in certain regions of magnetic configuration space. Furthermore, by 
averaging over constant energy trajectories, one can construct an equivalent dynamical flow 




(8) 




1 




6 



equation in energy space. This kind of approach has aheady been used in the hterature^ and 
can lead to the appearance of stable limit cycles at currents less than the critical current 
as can be intuitively inferred by considering which constant energy trajectories lead to a 
canceling of the flow in (11). 

Starting from an initially stable magnetic state, spin-torque effects will tend to drive the 
magnetization toward the current's polarization axis. Once the current is turned off, the 
projection of the magnetization vector along the uniaxial anisotropy axis will almost surely 
determine which stable energy state (parallel or anti-parallel) the magnetic system will relax 
to as long as the energy barrier ^ is large enough. As such, switching dynamics are best 
studied by projecting equations (8) along the uniaxial anisotropy axis uk- One then obtains 
a stochastic differential equation describing the dynamics of such a projection: 



q = a 



[nj + g)(l - q^) + njq{q -) 



+ aHn,my + ^j{l-q^)oW. (10) 

In the above equation, we have defined q = m ■ nx', and are the projections of 
the uniaxial anisotropy axis respectively on the z and x axes. Furthermore, the multiple 
stochastic contributions are assumed to be gaussian random variables with identical average 
and space dependent variance, ly is a standard mean zero, variance 1, Wiener process, and 
its prefactor explicitly expresses the strength of the stochastic contribution^. As will be 
shown in the following subsections, (12) is a convenient analytical tool in specific scenarios. 
In general, though, it is not useful as it explicitly depends on the dynamics of both the m^, 
and rriy components of the magnetization. 

Numerically, we can solve (8) directly. In all our simulations, the initial ensemble of 
magnetizations was taken to be Boltzmann-distributed along the anti-parallel orientation. 
We assume that the energy barrier height is so large that, before current effects are activated, 
thermalization has only been achieved within the antiparallel energy well and no states have 
had time to thermally switch to the parallel orientation on their own. A typical histogram 
of magnetic orientations at a given time is shown in Figure 1. 

We now turn on a current and allow the system to evolve for a fixed amount of time. 
Once this time has passed, we use the projection rule expressed above to evaluate what 
fraction of the ensemble has effectively switched from the anti-parallel to the parallel state. 
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FIG. 1: Histogram distribution of after letting the magnetic system relax to thermal equilibrium 
(10^ natural time units). The overlayed red dashed line is the theoretical equilibrium Boltzmann 
distribution. In the inset we show a semilog-plot of the probability vs. dependency. As 
expected, the data scales linearly. 

V. COLLINEAR SPIN-TORQUE MODEL 

Having derived the necessary expressions for our macrospin model dynamics, it is useful 
to consider the following simplification. Let us take the uniaxial anisotropy and spin-current 
axes to be coUinear, namely, uk = hp = k. In such a scenario, the stochastic LLG equation 
simplifies significantly. In particular, (11) reduces to the simplified form: 



a 



q = a{I + q){l - q^) + J^il - q^) oW. 



(11) 



In this symmetric scenario, q coincides with and magnetization reversal has been reduced 
to a straightforward 1-D problem. For a general value of / < 1, the evolution of q has two 
local minima and a saddle. The two stable configurations are at g = — 1 and q = 1, while 
the saddle is located at g = — /. For currents / > Jc = 1 there is only one stable minimum. 
Above critical current, spin torque pushes all magnetic configurations toward the parallel 
q = 1 state. 
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A. Collinear High Current Regime 



In the high current regime / ^ /c we expect the deterministic dynamics to dominate over 
thermal effects. We refer to this also as ballistic evolution interchangeably. The determistic 
contribution of (12) can then be solved analytically given an initial configuration q = rriz = 
—rriQ. The switching time Tg will simply be the time taken to get from some ruz = —mo < 
to rriz = and reads: 



1 /-o 
r,(mo) = - / 
a J-'i 



2a{P - 1) 



dm 


f m)(l 


— m^) 


{/log 


'1 + mo' 




-1 - niQ. 



log 



mr 



2 log 



(12) 



J -mo_ 

Since the magnetic states are considered to be in thermal equilibrium before the current is 
turned on, one should average the above result over the equilibrium Boltzmann distribution 
in the starting well to obtain the average switching time {ts)b- For ^ large enough, such an 
initial distribution will be: 



PB[m) 



exp[,^m^ 



(13) 



where F[x] = exp(— x^) exp{y'^)dy is Dawson's integral. This expression can be used to 
compute the average switching time numerically. 

As the intensity of spin-currents becomes closer to Jc, thermal effects increasingly con- 
tribute. Moreover, diffusion gradients add to the deterministic drift, which can be shown 
explicitly by writing (13) in its equivalent Ito form. Doing so leads us to a first correction 
of the ballistic dynamics due to thermal influences. The z-component behavior then reads: 



rhz = a{I + m^)(l — ml) 



(14) 



2e ■ V e 

The first term on the right hand side is still the ballistic flow that we have just discussed. 
The second term is the desired diffusion-gradient drift term. The contribution of such a term 
generates a net motion away from the stable minima of the ballistic equations as one expects 
to see under the influence of thermal effects. Again, we can solve the drift dominated flow 
analytically to compute the switching time. Considering diffusion-gradient drift, this reads: 
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FIG. 2: Blue line shows the fit of the ballistic limit to the numerical data (in blue crosses). Red 
line shows the improvement obtained by including diffusion gradient terms. Times are shown in 
units of (s ■ T) where T stands for Tesla: real time is obtained upon division by Hk- 



1 /-o 
rs(mo) = - / 

a J-r 



dm 

mo (/ + m)(l - m2) - (m/20 
1 X ^ \og[{wj —mo)/w 
a 



1^ iog[{Wj -mo)/Wji 

a^3w]-2Iw,-il-^) ^ ^ 

Where the wj are the three zeros of the cubic equation — Iw"^ — (1 — ^)'"^ + / = 0. As 
before, the average switching {ts)b time will simply be given by averaging numerically over 
the Boltzmann distribution pB- In Figure 2, the reader can see how well these two limiting 
results fit the simulation data. As expected both expressions coincide in the limit of high 
currents. 



VI. UNIAXIAL TILT 

In the high current regime (J ^ Jc), where nx = k (i.e. the uniaxial tilt is aligned with 
the z-axis), the ballistic equation for was shown to decouple from the other components, 
and the dynamics became one dimensional and deterministic. For the more general case 
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FIG. 3: q: green > 0, red < for applied current 1 = 5. The plane dissecting the sphere is 
perpendicular to the uniaxial anistropy axis. Its intersection with the sphere selects the regions 
with highest uniaxial anisotropy energy. 

where the uniaxial anisotropy axis may have any tilt with respect to the z-axis, such a 
critical current is not as intuitively defined. Unlike the coUinear limit, a critical current, 
above which all magnetic states perceive a net flow towards an increasing projection, does 
not exist. One can plot q = ih ■ Uk over the unit sphere to see what regions allow for an 
increasing and decreasing projection as the current is changed. An example of this is shown 
in Figure 3. 

Unfortunately, regions characterizing negative projection flow can be shown to persist 
at all currents. The approach is reflned by requiring that on average, over constant-energy 
precessional trajectories, the flow is toward the positive uniaxial anisotropy axis: (rii-n/^) > 
0. Such trajectories are found by solving the flow equations with = a = 0. Solutions 
correspond to circular libations about the uniaxial anisotropy axis. The critical current is 
then redeflned to be the minimum current at which the average projectional flow is positive 
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(t)(T.s) 

FIG. 4: Mean switching time behavior for various angular tilts above critical current. Each set 
of data is rescaled by its critical current such that all data plotted has /c = 1. Angular tilts are 
shown in the legend in units of 7r/36 such that the smallest angular tilt is and the largest is vr /4. 
Times are shown in units of {s ■ T) where T stands for Tesla: real time is obtained upon division 

at all possible precessional energies. This is easily done and results in: 



/ > max. 



— e 



cos u 



1 



COS w 



'■crity 



(16) 



thus allowing for a direct comparison of dynamical switching results between different angu- 
lar configurations of uniaxial tilt. In our discussion of (12) we mentioned how in the general 
case, presenting uniaxial tilt, there is no way to reduce the dimensionality of the full dy- 
namical equations. In fact, in the presence of tilt, precessional trajectories might allow for a 
magnetization state to temporarily transit through a g > 0, "switched" configuration, even 
though it might spend the majority of its orbit in a g < 0, "unswitched" configuration. This 
allows for a much richer mean switching time behavior, especially above critical current, as 
shown in Figure 4 and discussed more in depth later. 
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VII. THERMALLY ACTIVATED REGIME 

For currents I < Ic the switching rehes on thermal effects to stochastically push the 
magnetization from one energy minimum to the other. It is of interest to understand how 
switching probabilities and switching times depend on temperature and applied current. 
This is easily done in stochastic systems with gradient flow. In such cases an energy land- 
scape exists and a steady state probability distribution can be constructed to be used via 
Kramer's theory in deriving approximate low-noise switcing probabilities. 

Unfortunately, though, spin-torque effects introduce a non-gradient term, and the result- 
ing LLG equation does not admit an energy landscape in the presence of applied current. 
The collinear simplification, however, is an exception. As already described, in the absence 
of uniaxial tilt the dynamics become effectively one dimensional since the component 
decouples from the other magnetization components. Consider then (13): because it is de- 
coupled from the other degrees of freedom, we can construct a corresponding one-component 
Focker-Planck equation. The evolution in time of the distribution of q is then: 



For high energy barriers and low currents, the switching events from one basin to the other 
are expected to be rare. The probability of a double reversal should be even smaller. We 
therefore model the magnetization reversal as a mean first passage time (MFPT) problem 
with absorbing boundaries at the saddle point. The MFPT will then be given by— the 



dtpiq,t) = L[p\{q,t), 



(17) 



where 




(18) 



solution of the adjoint equation {L^T){q) 



1): 




(19) 



subject to the boundary condition (t)(0) 



0. This can be solved to give: 




(20) 



13 




Data 

Exponent = 1 
Exponent = 2 



0.85 



0.30 



0.70 



0.65 



10 



10 



10 



10 



(t)(T-s) 



FIG. 5: Mean switching time behavior in the sub-critical low current regime. Times are shown 
in units of {s ■ T) where T stands for Tesla: real time is obtained upon division by Hk- The red 
and green line are born by fitting to the data the functional form (r) = C exp(— ^(1 — I)^)-, where 
/i is the debated exponent (either 1 or 2) and C is deduced numerically. The red curve fits the 
numerical data asymptotically better the green curve. 

The rightmost integral can be computed explicitely. Retaining only dominant terms, the 
final integral can be computed by saddlepoint approximation to give: 



Such a square exponential dependence has recently been derived by Taniguchi and Ima- 
murap^i^, although a r oc exp(,^(l— /)) dependence, proposed elsewhere in the literature^i^ii^, 
has also been successfully used to fit experimental datai^. 

To decide between these experimental dependences, we fit the scaling behaviors in Figure 
5, along with the theoretical prediction from (24). The square exponential dependence 
fits the data better, confirming our analytical results. Furthermore, comparison to the full 
theoretical prediction demonstrates that even for mean switching times of the order 10~^ 
milliseconds, asymptoticity still is not fully achieved. 

All that remains is to consider the effects of angular tilt on the switching properties in 




v/?exp(e(l-/mv/e(l-/)) 



(21) 
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the thermally activated regime. Insight into this problem can be obtained by invoking (12) 
again. For small values of a, the term in square brackets is of leading order over the second 
ballistic term depending on rriy. This allows us, in the small a regime, to neglect the second 
ballistic term altogether. 

We now concentrate on the behavior of the term in square brackets. For low sub-critial 
currents, switching will depend on thermal activation for the most part. We expect an 
initially anti-parallel configuration to not diffuse very far away from its local energy minima. 
We expect it to remain that way until a strong enough thermal kick manages to drive it across 
the energy barrier. Because of this, the second term appearing in the square brackets will 
generally be close to zero as the particle awaits thermal switching. To make the statement 
more precise, one can imagine the magnetic state precessing many times before actually 
making it over the saddle. The second term can then be averaged over a constant energy 
trajectory and shown to vanish identically. Hence, in the subcritical regime, (12) can be 
rewritten in the following approximate form: 



This, is reminiscent of the ID projectional dynamics discussed in relation to the collinear 
limit, and shown explicitly in (13). The only difference between the two is the substitution 
/ — !■ rizl; recall because = cos{uj), Uzl = I/Ic- In other words, the thermally activated 
dynamics are the same for all angular separations up to a rescaling by the critical current. 
We then expect that the mean switching time dependences remain functionally identical to 
the collinear case for all uniaxial tilts. We have confirmed this by comparison with data from 
our simulations, and the results are shown in Figure 6. As predicted, all mean switching 
time data from different uniaxial tilts collapses on top of each other after a rescaling by each 
tilt's proper critical current. 

In comparing our scaling relationships between current and mean switching time with 
the previous literature, a subtle issue must be addressed. Results obtained by Apalkov and 
Visscher- rely upon an initial averaging of the dynamics in energy space over constant energy 
trajectories (limit for small damping) and only subsequently applying weak noise methods 
to extrapolate switching time dependences. The small damping and weak noise limits are 
singular and the order in which they are taken is important. Both limits radically alter the 
form of the equations: whereas both limits suppress thermal effects, the first also severely 




(22) 
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FIG. 6: Mean switching time behavior in the sub-critical low current regime. Various uniaxial tilts 
are compared by rescaling all data by the appropriate critical current value. Times are shown in 
units of (s • T) where T stands for Tesla: real time is obtained upon division by Hk- 

restrains the deterministic evolution of the magnetic system. Our approach considers the 
weak noise Umit and, only in discussing the effects of an angular tilt between polarization 
and easy axes do we employ the small damping averaging technique to obtain functional 
forms for the mean switching time. The switching time data shown seems to justify, in this 
particular, an interchangeability between these two limits. More generally, however, one 
should not expect the two limits to commute. 

VIII. SWITCHING TIME PROBABILITY CURVES 

Up until now, we have analyzed the main properties of spin-torque induced switching 
dynamics by concentrating solely on the mean switching times. In experiments, though, one 
generally constructs full probability curves. The probability that a given magnetic particle 
has a switching time Tg < t can be explicitely written as: 
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/••m{T) 

P[rs<T]= / dxpsix 
Jo 



= expK(l - mirY)] '^^^ '\ (23) 

where m(r) is the initial magnitization that is switched deterministically in time r. Once 
one has evaluated m(r), the probability curve follows. Ideally, in the ballistic regime, one 
would like to invert the ballistic equations. Unfortunately, though, the solutions of such 
ballistic equations are generally transcendental and cannot be inverted analytically. Even 
in the simpler coUinear case, as can be seen from equations (14) and (18), no analytical 
inversion is possible. One must instead compute the inversion numerically—. Nonetheless, 
one can construct appropriate analytical approximations by inverting the dominant terms 
in the expressions. In the case of (14), for example, one has that for currents much larger 
than the critical current: 

t \ -^1 i^ + rn . . 



which can be inverted to give: 



P - 1 

m(r) = tanh[ar — j — ]. (25) 



Plugging into expression (25) for the r probability curve, one has: 

F[v/etanh[ar^]] 



-P[ts < t] = exp 



cosh^ [ar^^^^ 

This expression can be truncated to a simpler form by noting that the leading exponential 
term dominates over the ratio of Dawson functions. Furthermore, if one considers the limit 
of large values for r (or, analogously, / ^ 1), the 'cosh' can be also approximated by its 
leading exponential term. We are finally left with: 



(26) 



P[ts < t] — exp 



72-11 



cosh [«r— 



exp 



-4^ exp[— 2ar- 



(27) 
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which is very similar in form to what has aheady been reported and used for fitting in the 
hterature^ii^. 

In the low current regime, one constructs probability curves by considering the mean 
switching time and modeling a purely thermal reversal as a decay process with rate given 
by equation (24). The fraction of switched states then vary in time as: 



Upon introducing uniaxial tilt, precessional effects can be witnessed directly on the switching 
probability curves in the super-critical regime. One expects that in the initial phases of 
switching, the fraction of switched states is sensitive to the time at which the current is 
turned off. One may accidentally turn off the current during a moment of transient passage 
through the switched region along the precessional orbit. This was checked and verified from 
our numerical simulations (see Figure 7). More generally, effects similar to the "waviness" 
seen in the mean switching time curves can be seen in the probability curves as well as the 
angle of uniaxial tilt is allowed to vary (see Figure 8). Only a numerical solution of the LLG 
equation can bring such subtleties to light. 

IX. CONCLUSION 

We have constructed the theory underlying the dynamics of a uniaxial macrospin subject 
to both thermal fluctuation and spin-torque effects. We then studied the subtle interplay 
between these two effects in aiding magnetization reversal between energy minima in a 
magnetically bistable system. Two regimes stand out in such a theory: a ballistic regime 
dominated by the deterministic flow and a thermally activated regime where reversal is 
dominated by noise. In the ballistic regime we discussed how to approximate the mean 
switching time behavior and found that corrections due to the diffusion-gradient term, arising 
from the stochastic equations, allow one to model the dynamics more accurately. 

In the thermally activated regime, we solved the relevant mean first passage time problem 
and obtained an expression for the dependence of mean switching time on applied current. 
In doing so, the correct scaling was shown to be (r) oc exp(— ,^(1 — /)^), in contrast to 
the prevailing view that (r) oc exp(— ^(1 — /)). Analytical results were then compared 




(28) 
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Switching Probability as a Function of Time 




25 30 35 40 45 50 
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FIG. 7: Influence of precessional orbits on transient switching as seen from the switching time 
probability curve in the supercritical current regime. The case shown is that of an angular tilt of 
7r/4 subject to a current intensity of 1.5 times the critical current. The non-monotonicity in the 
probability curve shows the existence of transient switching. 

with detailed numerical simulations of the stochastic LLG equation. The massively parallel 
capabilities of our GPU devices has allowed us to explore the behavior of macrospin dynamics 
over six orders of temporal magnitude. Comparing to our analytical results, we suggest 
that the thermal asymptotic behavior is achieved very slowly with respect to the switching 
timescales that have been probed experimentally. 

Different geometrical configurations of the uniaxial anisotropy axis with respect to the 
spin-current axis were shown to influence the thermally activated regime very minimally 
inasmuch as the currents were rescaled by the proper critical current of the angular setup. 
Only in the super-critical regime were distinctions shown to exist due to complex precessional 
and transient switching behavior. 

These results have important implications for the analysis of experimental data in which 
measurements of the switching time versus current pulse amplitude are used to determine 
the energy barrier to magnetization reversal. Clearly use of the correct asymptotic scaling 
form is essential to properly determine the energy barrier to reversal. The energy barrier, 
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Time 



FIG. 8: Spin-torque induced switching time probability curves for various angular configurations 
of uniaxial tilt (a sample normalized current of 10 was used). 

in turn, is very important in assessing the thermal stability of magnetic states of thin film 
elements that are being developed for long term data storage in STT-MRAM. Further work 
should address how these results extend to systems with easy plane anisotropy and situations 
in which the nanomagnet has internal degrees of freedom, leading to a break down of the 
macrospin approximation. 

We also note that current flow is a source of shot noise, which at low frequencies acts like 
a white-noise source in much the same way as thermal noise. It is therefore interesting to 
understand when this additional source of noise plays a role. For a magnetic layer coupled 
to unpolarized leads, the current induced noise on the magnetization dynamics was found 
to be j^^^^jj^^V—, where V is the voltage drop across the magnetic layer, while Fi/F/j is a 
dimensionless ratio characterizing the coupling strength of the magnetic layer to the left (L) 
and right (R) leads. Thus the noise is maximal (V/A) for perfectly symmetrical couplings, 
and is smaller in the limit of highly asymmetric contacts. This basic behavior, and the order 
of magnitude of the effect, is not likely to be modified by polarized leads. We argue that the 
temperatures at which the experiments have been performed current noise effects are not 
important. The experiments are performed at room temperature where T = 300 K. For an 
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all-metallic device, such as a spin-valve nanopillar, the couplings are nearly symmetrical and 
at the critical current a typical voltage drop across the magnetic layer is less than 10 yuV or, 
equivalently, 1 K. For a magnetic tunnel junction device V can be ~ 1 V. However, in this 
case the coupling is asymmetric. One lead (L) forms a magnetic tunnel junction with the 
nanomagnet, while the other (R) a metallic contact. This gives F^/F/^ > 10^ and a relevant 
energy ~ 1 K, again far lower than room temperature. It appears that current induced 
noise can only be important at room temperature for a nanomagnet coupled symmetrically 
between two tunnel barriers. 
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